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Drawing from exact, approximate and numerical results an overview of the properties of the 
out of equilibrium response function in phase ordering kinetics is presented. Focusing on the zero 
^s^j . field cooled magnetization, emphasis is on those features of this quantity which display non trivial 

behavior when relaxation proceeds by coarsening. Prominent among these is the dimensionality 
t dependence of the scaling exponent a x which leads to failure of the connection between static and 

. dynamic properties at the lower dimensionality di, where a x = 0. We also analyse the mean 

spherical model as an explicit example of a stochastic unstable system, for which the connection 
0^ ' between statics and dynamics fails at all dimensionalities. 

<n : 

PACS numbers: 
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(1) ■ INTRODUCTION 

It is a great pleasure to present this talk on the occasion of the 60th birthday of Francesco Guerra. My interest 
in relaxation phenomena, in fact, initiated quite a few years ago in Salerno where Nelson stochastic mechanics was 
a subject very much cultivated under the guide of Francesco. I have benefited a lot from his teaching and from the 
^— 1 . very stimulating atmosphere he was able to create about the applications of stochastic methods to physics in general. 

Coming to the subject of this talk, phase ordering is usually regarded as the simplest instance of slow relaxation, 
quite useful for a first and easy to understand approach to concepts like scaling and aging, which are the hallmarks 
of glassy behavior 0. However, next to the similarities there are also fundamental differences which require to 
keep phase ordering well distinct from the out of equilibrium behavior in glassy systems, both disordered and non 
disordered. The main source of the differences is the simplicity of the free energy landscape in the case of phase 
ordering compared to the complexity underlying glassy behavior. This leads to a simple low temperature state for 
phase ordering systems as opposed to the complexity of replica symmetry breaking for glasses (at least in the mean 
field picture of spin glasses). An analogous sharp distinction is believed to exist also in the behavior of the out of 
equilibrium linear response function. As a matter of fact, in the framework where static and dynamic properties are 
connected systems may be classified on the basis of the fluctuation dissipation relation 0. 

In the jargon of slow relaxation, phase ordering is frequently referred to as coarsening. So, if coarsening is associated 
to simplicity and triviality of behavior what is there to investigate about it? Besides the obvious motivation that the 
basic, paradigmatic cases need to be thouroughly understood, an additional reason, among others, is that in some 
cases the existence of complex slow relaxation is identified through the exclusion of coarsening. This requires precise 
and reliable knowledge of all what goes on when relaxation proceeds by coarsening. An example comes from the long 
standing controversy about the nature of the low temperature phase of finite dimensional spin glasses. One recent 
argument in favour of replica symmetry breaking is that the observed behavior of the response function is incompatible 
with coarsening 0,0. This might well be the case. However, for the argument to be sound, the understanding of 
the out of equilibrium behavior of the response function during phase ordering needs to be up to the level that such 
a delicate issue demands. 

It is the purpose of this talk to present an overview of the accurate investigation of the response function in phase 
ordering that we have carried out in the last few years. Focusing on the zero field cooled magnetization (ZFC), which 
is one instance of integrated response function, it will be argued that the response function in phase ordering systems 
is not as trivial as it is believed to be and, after all, it is not the quantity best suited to highlight the differences 
between systems with and without replica symmetry breaking. In fact, as we shall see, phase ordering and therefore 
a replica symmetric low temperature state are compatible with a non trivial ZFC. When this happens there is no 
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connection between static and dynamic properties. Phase ordering systems offer examples of two distinct mechanism 
for the lack of this important feature of slow relaxing systems, the vanishing of the scaling exponent of ZFC and 
stochastic instability. 



PHASE ORDERING 



Let us first briefly recall the main features of a phase ordering process. Consider a system, like a ferromagnet, 
with order parameter (vector or scalar, continous or discrete) cj>(x) and Hamiltonian TL[<j){x)\ such that below the 
critical temperature Tq the structure of the equilibrium state is simple. For example, in the scalar case, there are 
two pure ordered states connected by inversion symmetry. The form of the Hamiltonian can be taken the simplest 
compatible with such a structure, like Ginzburg-Landau-Wilson (GLW) for continous spins or the nearest neighbors 
Ising Hamiltonian for discrete spins. 

In the following we will be interested in the space and time dependent correlation function 



C(r,M) - (mt)^,s)} - (4>{x,t)){4>{x',s)) 



(1) 



where the average is taken over initial condition and thermal noise, r = x — x' and t > s > are two times after 
the quench. The conjugated linear response function is given by 



R(f,t,s) 



5(<f>(x,t)) 



Sh(x', s) 



(2) 



h=0 



and ZFC is defined by 



X(r,t,t w ) = J dsR(r,t,s). 



(3) 



Dynamics over phase space: equilibration versus falling out of equilibrium 

For a temperature T below Tc, in the thermodynamic limit, the phase space £1 = {[(f>(x)]} may be regarded as 
the union of three ergodic components = il + U £1_ U fto, where Q± and fig are the subsets of configurations 
with magnetization limy^oo X fy dx(f>{x) positive, negative and vanishing, respectively. Denoting by p± [<t>(x)\ the 
two broken symmetry pure states, all equilibrium states are the convex linear combinations of p±. In particular, the 
Gibbs state is the symmetric mixture pg[4>(x)] — J? exp(— Tt[4>(x)]/T) = ^p + [(j)(x)) + ^p-[(j>(x)]. The Q± components 
are the domains of attraction of the pure states with p + (fl + ) = = 1 and f^o is the border in between them, 

with zero measure in any of the equilibrium states. 

When ergodicity is broken, quite different behaviors may arise depending on the initial condition po[cj)(x)] = 
p([4>(x)],t = 0). Here, we consider the three cases relevant for what follows, assuming that there are not explicit 
symmetry breaking terms in the equation of motion: 

1. equilibration to a pure state 

if po(Q+) = 1 or po(Q-) = 1, in the time evolution configurations are sampled from either one of f2± and 
p([(j)(x)],t) equilibrates to the time independent pure state p±[4>{x)] within the finite relaxation time t eq ~ £ z , 
where £ is the equilibrium correlation lenght and z is the dynamic exponent entering the domain growth law, to 
be defined shortly. The correlation function is the same in the two ergodic components and, after equilibration, 
is time translation invariant 



C st (r,t- s) = (cf>(x, t)4>(&, s)) ± - M 2 (4) 

where {4>{x))± — ±M is the spontaneous magnetization. For large distances r > ( an< i time separations 
t — s ^> t eq , the clustering property {4>(x, t)(f)(x' , s)}± — {4>(x, t))±(4>(x' , s)}± is obeyed and the correlations 
decay to zero, as required by ergodicity. 
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2. equilibration to the Gibbs state 

if po(£l+) = po(£l-) = 1/2, then configurations are sampled evenly from both disjoint components f2+ and f2_. 
The probability density p([(f>(x)],t) equilibrates now to the Gibbs state pg[4>{x)\ with the same relaxation time 
t eq as in the relaxation to the pure states. Broken ergodicity shows up in the large distance and in the large 
time properties of the correlation function. After equilibration, one has 

C G (f, t-s) = C st (r, t-s) + M 2 (5) 
from which follows that correlations do not vanish asymptotically or that the clustering property is not obeyed 

lim C G (r,t-s)= lim C G (r,t-s) = M 2 . (6) 

r— >oo (t-s)^oo 

3. falling out of equilibrium over the border 0,0/ 

if po(^o) = 1> f° r the infinite system p(f2o,i) = 1 also at any finite time after the quench. Namely, the system 
does not equilibrate since in any equilibrium state the measure of f2o vanishes. Phase ordering corresponds 
to this case. In fact, the system is initially prepared in equilibrium at very high temperature (for simplicity 
Tj = oo) and at the time t = is suddenly quenched to a final temperature T below Tq- In the initial state 
the probability measure over phase space is uniform po[(f>(x)] — implying that the initial configuration at 

t = belongs almost certainly to f^o, since with a flat measure |Oq| is overwhelmingly larger than |f2±|. 

The morphology of typical configurations visited as the system moves over fio is a patchwork of domains of 
the two competing equilibrium phases, which coarsen as the time goes on. The typical size of domains grows 
with the power law L(t) ~ t 1 ^, where z = 2 (independent of dimensionality) for dynamics with non conserved 
order parameter 0, as it will be consider here. The sampling of configurations of this type is responsible of 
the peculiar features of phase ordering. At a given time t w there remains defined a length L(t w ) such that 
for space separations r <C L(t w ) or for time separations t — t w <C t w intra-domains properties are probed. 
Then, everything goes as in the case 2 of the equilibration to the Gibbs state, ergodicity looks broken and the 
correlation function obeys Eq. J5|. Conversely, for r 3> L(t w ) or t/t w 3> 1, inter-domains properties are probed, 
ergodicity is restored (as it should be, since evolution takes place within the single ergodic component f^o) an d 
eventually the correlation function decays to zero. However, the peculiarity is that if the limit t w — > oo is taken 
before r — * oo, in the space sector ergodicity remains broken giving rise, for instance, to the growth of the Bragg 
peak in the equal time structure factor. 

According to this picture, the correlation function can be written as the sum of two contributions 

C(r, t, s) = C st {f, t-s) + C ag (r, t, s) (7) 

where the first one is the stationary contribution of Eq. Q describing equilibrium fluctuations in the pure states 
and the second one contains all the out of equilibrium information. The latter one is the correlation function 
of interest in the theory of phase ordering where, in order to isolate it, zero temperature quenches are usually 
considered as a device to eliminate the stationary component. It is now well established that C ag {r, t, s) obeys 
scaling in the form [2j 

C ag (f,t,s)=d(r/L(s),t/s) (8) 
with C(x, y) = M 2 for x < 1 and y ~ 1, while 



d(r/L(s),t/s)~(t/s)-V°h(r/L(s)) 
for large time separation 0, where A is the Fisher-Huse exponent. 



(9) 
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ZERO FIELD COOLED MAGNETIZATION 



Let us next consider what happens when a time independent external field h(x) is switched on at the time t w . To 
linear order the expectation value of the order parameter at the time t is given by 

((f>(x,t)) h = ((f>(x,t))o+ [ dx' X (x-x',t,t w )h(x') (10) 



and if h(x) is random with expecations h(x) = 0, h(x)h(x') — ft.Q<5(x — x') then one has 



X(x - y,t,t w ) = h (cj)(x,t)}hh{ff). (11) 

Namely, ZFC is the correlation at the time t of the order parameter with the random external field. 
Going to the three processes considered above 

1. after equilibration in the pure state has occurred and the stationary regime has been entered, the order parameter 
correlates with the external field via the equilibrium thermal fluctuations. The fluctuation dissipation theorem 
(FDT) is obeyed 

Xst{r,t-t w ) = I[C at (r,t- s = 0)-C st (r,t-t w )} (12) 
and since C st (r, t — t w ) decays to zero for t — t w > t eq , over the same time scale Xst(r, t — t w ) saturates to 

lim Xst(r,t-t w ) =Xeq{r) = ^C eq {f) (13) 

t — >oo ± 

which is the susceptibility computed in the final equilibrium state p±[(j}{x)]. 

2. As far as ZFC is concerned, there is no difference between the relaxation to the mixed Gibbs state and the 
relaxation to a pure state. Hence, FDT is satisfied and can be written both in terms of C s t or Cg since, as 
Eq. JSJ) shows, they differ by a constant. 

3. In the phase ordering process the system stays out of equilibrium, so it useful to write ZFC as the sum of two 
contributions |9|| 



X(r, t,t w ) = Xst(r,t- t w ) + Xag(r 7 t,t w ) (14) 

where Xs*( 7? j t ~ tw) satisfies Eq. (|12f) and Xag(^ t, t w ) represents the additional out of equilibrium response. In 
connection with this latter contribution there are two basic questions 

i) how does it behave with time 

ii) what is the relation between Xag and C agi if any. 



Scaling hypothesis 



Since ZFC measures the growth of correlation between the order parameter and the external field, the first question 
raised above addresses the problem of an out of equilibrium mechanism for this correlation, in addition to the thermal 
fluctuations accounting for Xst- Restricting attention from now on, for simplicity, to the case of coincident points 
(r = 0) and dropping therefore the space dependence, the starting point for the answer is the assumption of a scaling 
form 
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Xag(t,t w ) ~ t- a -X ag {t/t w ) (15) 

which is the counterpart of Eq. © for the correlation function. 

The next step is to make statements on the exponent a x and on the scaling function X ag (x) ■ There exists in the 
literature an estimate of a x based on simple reasoning. What makes phase ordering different from relaxation in the 
pure or in the Gibbs state is the existence of defects. The simplest assumption is that Xag{t,t w ) is proportional to 
the density of defects 0, 0, . This implies 

a x = 5 (16) 
where the exponent 8 regulates the time dependence of the density of defects Pdefect(t) ~ L(t)~ n ~ t~ 5 , namely 



5 = n/z (17) 

with n = 1 for scalar and n = 2 for vector order parameter [l|. 

According to this argument a x should be independent of dimensionality. This conclusion is not corroborated by 
the available exact, approximate and numerical results. On the basis of exact analytical solutions for the d = 1 Ising 
model 0, 0| and for the large N model 0] , approximate analytical results based on the gau ssian auxiliary field 
(GAF) approximation [13, and numerical results from simulations of the Ising model [ill [ili Il7j | with d = 2,3, 4, 
the following general formula for a x has been obtained 



_j 6 {jA) fOT d<d " 

a x ~ \ 6 with log corrections for d — du 0-°) 
[ 6 for d> du 

where du and du > du do depend on the system in the following way 

• du is the dimensionality where a x = 0. In the Ising model du — 1, while in the large N model du = 2. The 
speculation is that in general du — 1 for systems with discrete symmetry and du = 2 for systems with continous 
symmetry, therefore suggesting that du coincides with the lower critical dimensionality of equilibrium critical 
phenomena, although the reasons for this identification are far from clear. 

• du is a value of the dimensionality specific of ZFC and separating d < du, where a x depends on d, from d > du 
where a x is independent of dimensionality and Eq. I|16fl holds. The existence of du is due 01 to a mechanism, 
i.e. the existence of a dangerous irrelevant variable, quite similar (including logarithmic corrections) to the one 
leading to the breaking of hyperscaling above the upper critical dimensionality in static critical phenomena. 
However, du cannot be identified with the upper critical dimensionality since we have found, so far, du = 3 in 
the Ising model and du = 4 in the large N model. In the scalar case it may be argued 0] that du coincides 
with the dimensionality du = 3 such that interfaces do roughen for d < du and do not for d > d^. A general 
criterion for establishing the value of du is not yet known. 

The validity of Eq. (|15fl with a x given by Eq. I |18| l has been checked, in addition to the cases where analytical results 
are available, with very good accuracy in the simulations of the Ising model 0|. The values of 8, du and du obtained 
for the different systems are collected in Table H] and the behavior of a x as dimensionality is varied is displayed in 

Figm 
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TABLE I: Parameters entering Eq. llBH in various models. 
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Eq.(18) 
Numerical 
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FIG. 1: Exponent a x in the Ising model at various dimensionalities. The continous line represents Eq. I|18^ . while the dots are 
the values from the exact solution [Tl|j of the model at d = 1 and from simulations at d = 2, 3, 4. 



At this point, a comment is in order. From Eqs. (|14fl and (|15fl follows that the existence of out of equilibrium degrees 
of freedom, or defects, generates the aging contribution to ZFC. However, when a x > this is bound to disappear 
as time becomes large, either t w — ► oo or t/t w — > oo. The vanishing of this contribution is fastest for d > djj, where 
a x = 5, and becomes slower and slower below du, as dh is approached. Nonetheless, as long as a x > 0, eventually 
Xag disappears and Eq. ijT3|) holds also for phase ordering processes for all d > d^. This is no longer true for d < d^ 
where a x < and, consequently, \ ag gives a contribution to the response which persists also in the asymptotic time 
region leading to 

lim X{t,tw) > Xeq (19) 
t— >oo 

where \ eq stands for Xeq(r = 0). In the following we will consider d> dL- 

Fluctuation dissipation relation 

Let us now come to the second question. A very useful tool for the study of slow relaxation phenomena has been 
introduced by Cugliandolo and Kurchan [19| through the off equilibrium fluctuation dissipation relation. This can 
be introduced in the following way. Given that C(t,t w ) is a monotonously decreasing function of t, for fixed t w it is 
possible to invert it and write 



x(t,t w ) = x(C(t,t w ),t w ). 

Then, if for a fixed value of C(t, t w ) there exists the limit 



(20) 



lim x(C,t w ) = S(C) 



(21) 



the function S(C) gives the fluctuation dissipation relation. In the particular case of equilibrium dynamics, FDT is 
recovered and 5(C) = [C(0) — C]/T. Originally introduced in the study of the low temperature phase of spin glas s 
mean- field models, the fluctuation dissipation relation has been found in many other instances of slow relaxation |20|. 

Now, in order to search for 5(C) in the case of phase ordering, let us set f = in Eq. (JHJ and let us eliminate t/t w 
between Xag and C ag obtaining 



Xag{tj tw) ~ t. u . x Xag{C a g) ■ 

Then, from Eqs. (|14I12I22H one can write the general relation 



(22) 
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X (t,t w ) = - [C st (0) - C st (t -t w )]+ t- a -x ag (C ag ). (23) 

Using the identity [C st (0) — C st {t — t w )] = [C st (0) + M 2 — C st (t — t w ) — M 2 ] and considering that in the time 
interval where C st {t~t w ) ^ 0, i.e. for short times, one can replace C ag (t/t w ) with M 2 or equivalently C st (t— t w )+M 2 = 
C(t,t w ), the above equation can be rewritten as 



X (t,t w ) = Xst(C) + t- a -Xag{C ag ) (24) 

where the function Xst(C) is defined by 

_ / [C(0) - C(t, t w )} for M 2 < C < C(0) , . 

TXst{C) -\[C(Q)-M 2 ] for C<M 2 . (25) 

Therefore, from Eq. (|24() we have, first of all, that for phase ordering systems the fluctuation dissipation relation 
exists only if a x > (i.e. for d > cIl) and that for a x > 



S(C)=x st (C). (26) 

For a x — a little more care is needed. Equation lt^|l yields x(Mt«) — Xst{C) + Xag(C ag )- Recalling that a x = 
occurs at d = d^, which coincides with the lower critical dimensionality, in order to have a phase ordering process a 
quench to T = is required. This, in turn, implies C s t(t, s) — and C ag (t, s) — C(t, s). Therefore, using Eq. Ij25(l we 
have 



S(C)= X : q +Xa 9 (C) (27) 

where xt q — h m T^o[C(0) — M 2 ]/T is the T = equilibrium susceptibility, which vanishes for hard spins while is 
different from zero for soft spins. Physical implications of these results are discussed in the next section. 



STATICS FROM DYNAMICS 



One of the main reasons of interest in the fluctuation dissipation relation is that it may provide a link between 
static and dynamic properties. This was first found by Cugliandolo and Kurchan for mean-field spin glasses and 
then established in general by Franz et al. for slowly relaxing systems. 

Let us first introduce the overlap probability function P{q) in the equilibrium state obtained when the perturbation 
responsible of x(i, t w ) is switched off. The question is how is P(q) related to the unperturbed overlap function P(q). If 
P(q) — P(q) the system is stochastically stable |2l|. A milder statement of stochastic stability is that P(q) coincides 
with P(q) up to the effects of a global symmetry which might be removed by the perturbation. In the Ising case, 
where the perturbation breaks the up-down symmetry, defining 

P{q) = 29{q)P{q) (28) 

the system is stochastically stable in the sense that P(q) = P{q)- If the system is not stochastically stable, P(q) is 
not related neither to P(q) nor to P(q). As we shall see, this is the case of the mean spherical model. 
The statement is that 

1. if S(C) exists 

2. if lim t _,oo x(t, tw) = Xeg 



3. if the system is stochastically stable 
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then one has 



-T 



d 2 S{C) 



dC 2 



= P(q)- 



(29) 



C=q 



We may now check if this chain of connections applies to phase ordering. In replica symmetric low temperature 
states, as for instance in ferromagnetic systems, the overlap function is always trivial and we have 



P(q) = l[5(q~M 2 )+5(q + M 2 )} 



(30) 



with 



P(q) = P(q)=6(q~M 2 ). (31) 
Computing the derivative in the left hand side of Eq. (|29[) and using Eqs. (|26[l and l|27|) . for d > we find 



-T 



d 2 S{C) 



dC 2 



5(q - M 2 ). 



(32) 



C=q 



Therefore, Eq. (|2*9"|) is satisfied. 



Failure at d = dh 



Let us, next, go to the case d = d^. For the sake of definiteness, we consider the case of the Ising model with 
d = 1. In order to make compatible the two requirements of having an ordered equilibrium state and a well defined 
linear response function, instead^of taking the T — > limit it is necessary to take the limit of an infinite ferromagnetic 
coupling PJ. Then, P(q) and^(g) are given by Eqs. and l|31|l with M 2 — 1 at all temperatures. On the other 



hand, for any T we also have |l lj 



TXag(C) 



V2 



arctan 



a/2 cot (| 



G 



(33) 



This gives 



-T 



d 2 S(C) 



dC 2 



C=q 



7rcos(7rg/2) sin(7rg/2) 
[2 - sin(7rg/2)] 2 ' 



(34) 



Hence, it is clear that Eq. I|29|l is not verified. The reason is that the second of the above conditions required for 
establishing the connection is not satisfied. In fact, from Eqs. I|24(l and (|33|l . keeping in mind that the limits t — > oo 
and C — > are equivalent, we have 



lim T X (t,t w ) = 1/V2 



(35) 



which is responsible of Eq. (|19f) and, therefore, of the violation of condition (2) above Eq. I|29|l . since in this case 

Xeq = 0. 



Failure by stochastic instability 

An interesting example [23j | , where statics cannot be reconstructed from dynamics because the third requirement of 
stochastic stability is not satisfied, comes from the spherical model. More precisely, one must consider in parallel the 
original version of the spherical model (SM) of Berlin and Kac [24| and the mean spherical model (MSM) introduced 
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by Lewis and Wannier |25|. with the spherical constraint treated in the mean. These two models are equivalent above 
but not below Tq |2rjj ]. The low temperature states are quite different, with a bimodal order parameter probability 
distribution in the SM case and a gaussian distribution centered in the origin in the MSM case. The corresponding 
overlap functions are also very different 23]. Considering, for simplicity, T — one has 




where Kq is a Bessel function of imaginary argument (Fig|3J). 




-10 1 q 



FIG. 2: Overlap distribution for mean spherical model with M 2 — 1. The arrows represent the 8 functions of the overlap 
distribution for the spherical model. 

However, after switching off an external held, one finds for both models P(q) = S(q — M ). This means that 
stochastic stability holds for SM but not for MSM. 

On the other hand, the relaxation properties are the same in the two models, both above and below Tc if the 
thermodynamic limit is taken before the t — ■+ oo limit 23). Then, the linear response function is the same for both 
models and obeys Eq. I)24|l with a x given by Eq. I|18|l . where 6, d^ and djj are the same as for the large N model 
(Table I). Hence, we have that although Eq. (|29J) is satisfied for both models, nonetheless statics and dynamics are 
connected only in the SM case, where P(q) = P{q). Instead, this is not possible in the MSM case where P(q) ^ P(q)- 

CONCLUSIONS 

We have shown that although the low temperature equilibrium state in phase ordering systems is simple, ZFC in 
quenches below Tc displays non trivial features. Summarising, these are 

• The aging component Xag{t,tw) obeys the scaling form 1|15|) with an exponent a x dependent on dimensionality 
according to Eq. (fT%|l . 

• The fluctuation dissipation relation S{C) is trivial, in the sense that it is consistent with a replica symmetric 
equilibrium state, for d > d^. However, as d — > d^ from above the non trivial contribution due to Xag(t,t w ) 
persists for longer and longer times as a x — > 0. 

• For d = dh the fluctuation dissipation relation is non trivial, since the aging contribution does not disappear 
asymptotically. From this follows i) that there is no connection between static and dynamic properties and 
ii) that on the basis of the behavior of ZFC is not possible to establish that the equilibrium state is replica 
symmetric. 

• The mean spherical model provides an explicit example of a system whose low temperature equilibrium state 
is not stochastically stable, producing the failure of the link between static and dynamic properties for all 
dimensionalities. 

• Extensive numerical simulations for ferromagnetic systems with scalar, vector, conserved and non conserved 
order parameter, at various space dimensionalities are under way. Preliminary results do show that the behavior 
of the response function above illustrated might well be generic for phase ordering systems |18| . 
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